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We investigate a t-J ladder model by numerical diagonalization 
method. By calculating correlation functions and assuming the Lut- 
tinger liquid relation, we obtained a possible phase diagram of the 
ground state as a function of J/t and electron density re. We also 
found that behavior of correlation functions seems to consist with the 
prediction of Luttinger liquid relation. The result suggests that the su- 
perconducting phase appear in the region of J/i>0.5 for high electron 
density and J/i>2.0 for low electron density. 
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Low- dimensional antiferomagnetic systems attract much interest due to 
the possible relevance to high-T c superconductivity. In particular, ladder 
systems which show a finite spin gap are extensively studied now. For exam- 
ple, recent measurements of the magnetic susceptibility and/or the nuclear 
relaxation rate of SrCu203 and (VO)2P2C>7 indicate that these materials are 
described weakly coupled ladder models with spin-gap |], 0, [|. Many theo- 
retical works also have been performed on Heisenberg of ladder models and 
the spin gap up to order of exchange energy is found @, ||, ||. Recently, It 
is claimed that the spin-gap survives against light hole doping and it van- 
ishes for heavy doping in t-J ladder models by numerical diagonalization 
method. 0, |[ |[ Hayward et al. gave numerical evidence of superconductiv- 
ity in t-J ladders for sufficiently large exchange energy J by a density matrix 
renormalization method |10|. 

If the low energy behavior of the ladder system can be described as that 
of one- dimensional (ID) system, an approach of the 'Luttinger liquid' theory 
will be useful. In the Luttinger liquid theory, some relations have been es- 
tablished as universal relations [11], [L2|. They supplied us with unambiguous 
and important information about not only one-band electronic systems but 



complicated 2-band systems [fl"3 



Hayward and Poilblanc |14| discussed the Luttinger liquid behavior of a 
t-J ladder model and give a possible phase diagram of a critical exponent 
K p which characterizes the algebraic decay of correlation functions in the 
Luttinger liquid theory. However, it is not clear whether or not exponents of 
correlation functions given by the Luttinger liquid relation consist with the 
behavior of correlation functions calculated by the numerical method [|lO| . 

In this letter, we reexamine the critical exponent and Luttinger liquid 
behavior of t-J ladder by numerical diagonalization method. We find that 
a prediction of the Luttinger liquid theory consists with the behavior of the 
correlation functions of t-J ladder. Using the exponent, we determine a 
superconducting region on the possible phase diagram. 

The t-J ladder Hamiltonian is written as 

H = -*£(4,a,<T C i+l,a,<T + h - c -) ~t 2(c!,i,«T c iAff + h - c -) 

i i 

+ J XX S '," ■ S i+llQ - ^n ii0[ n i)a ) + J " S *,2 - ^i,l n i,2) , (!) 

i,a i 

where c irT is the electron creation operator with spin a on a rung i of ladder 
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and Sj iQ is the spin operator made of c\ aa and q i(TiQ . Here, a(= 1,2) labels 

the two legs of the ladder and n i)Q , = J2a c l.a,a c i,a,a- J is exchange energy 
between the nearest neighbor sites and t stands for the transfer energy, which 
will be set to be unity hereafter in the present study. We take account of the 
infinite on-site repulsion by removing states with doubly occupied sites from 
the Hilbert space. 

We numerically diagonalize the Hamiltonian up to 20 sites (10 unit cells) 
using the standard Lanczos algorithm. We use the periodic boundary con- 
dition for N e = 4m + 2 and anti periodic boundary condition for N e = 4m, 
where N e is the total number of electrons and m is an integer. We also use 
the open boundary condition to calculate correlation functions because we 
are able to take more points in correlation functions. The filling n is defined 
by n = N e / (2N U ), where N u is the total number of unit cells ( each unit cell 
corresponds to a rung ). 

The chemical potential fi(N e , N u ) is defined by 

where E (N e ,N u ) is the total ground state energy. When the charge gap 
vanishes in the thermodynamic limit, the uniform charge susceptibility Xc is 
obtained from 

* (N « N - ) = W. + l,K)-l(N.-l,N.) ■ (3) 

It is noted that our definition of the \c is larger than that of previous works [§, 
nj by a factor of 2 since we use the charge susceptibility Xc per rung rather 



than per site. 

In the model which is isotropic in spin space, exponents of various types of 
correlation functions are determined by the critical exponent K p . For single 
band model, the Luttinger liquid theory predicts that Superconducting(SC) 
correlation function is dominant for K p > 1 (attractive case), whereas Charge 
or Spin Density Wave (CDW or SDW) correlation is dominant for K p < 1 
(repulsive case). In the case of non- interacting fermion systems, the exponent 
K p is always unity. 

On the other hand, the bosonization method |T5|, |T6| shows that two-chain 
model with a small interchain hopping has a spin-gaped phase and SC and 
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"4fcp" CDW correlations decay as ~ r 2K p and ~ r - 2K p respectively (SDW 
and "2/cf" CDW correlations decay exponentially). Hence, SC correlation is 
dominant for K p > 0.5. 

When we apply the above relations to the t-J ladder, we should pay 
attention to the filling of electrons in the t-J ladder. For low density region, 
we should adopt the relation for the single band model since electrons fill in 
only the lower band of the t-J ladder. For high density region, we can expect 
the relation for two-chain model to be useful. 

It is convenient to introduce a critical exponent K p which is related to 
the usual critical exponent K p by the relations , 

~ 7T 

K p = 2 VcXc ' ( 4 ) 

Vc = ^L( El -E ) , (5) 

where E\ — Eq is the lowest charge excitation energy. For low density, our 
definition of the K p is equivalent to that of the exponent K p fill , [I7fl . When 
electrons are filled in only lower band, we should define the K p of a multi- 



band model as above. Another definition of K p in previous works || [14 
is misleading in applying the Luttinger liquid relation to t-J ladder at low 
density. On the contrary, for high density, the K p is larger than the K p by a 
factor of 2 since we use the Xc given by eq.(3). Considering the \ c and K p 
of a non-interacting ladder system, it is easily understood. 

In both the low density and the high density cases, we expect that the 
superconducting correlations dominate for K p > 1. 

We can also determine the K p by the Drude weight D 

K P = n(^)^ , (6) 

= 

where E(<f>) is the total energy of the ground state as a function of flux 
0. Using these two independent equations of the K p , we can check the 
consistency of the Luttinger liquid relations. [0 

To investigate the band structure of the t-J ladder model, we show the 
chemical potential fi ( Fermi energy Ep at the T = ) as a function of electron 
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filling n in Fig.l. For a non- interacting ladder model, band structure is given 
as 

E ± (k) = 2tcos(k) ± t (8) 

where k is a wave vector. When n is smaller than 0.5, electrons are filled in 
only lower band E~(k). At n = 0.5, lower and upper bands begin to fill with 
electrons simultaneously. The slope of // is equal to zero and the value of the 
charge susceptibility Xc diverges at n = 0.5, where the slope of \i corresponds 
to 2/xc- Thus, the divergence of \ c signifies a change in the electronic state 
of non-interacting systems. 

On the other hand, Xc of t-J ladders seems to be finite at n — 0.5 for 
J<1.0. If the divergent point of Xc exists, it may move to more high density. 
In contrast to non-interacting fermion ladder, we expect that the ground 
state of t-J ladder continues from low density over n = 0.5. In fact, points of 
n = 0.5 for J<1.0 belongs to the spin-gapless region, which is yield by Ref. 
[§]. Roughly speaking, Ref.|§ shows that the spin-gap phase appears at high 
density and large J and the remainder is the gapless phase. 

It is interesting to see that the n— dependence of fi at J/t = 2 almost 
corresponds to that of non-interacting ladder system in low electron density. 
It suggests that the electronic state of non-interacting system resembles that 
of t-J ladder in this parameter region. To confirm the above, we calculate the 
density-density correlation function on the same leg Ccdw(R) =< ^o,a^i?,a > 
— < n 0tCt >< n# )Q > of both systems, where R is distance from one of ends 
of the system. In Fig. 2(a), we show Ccdw(R) for 18-site systems with 4- 
electrons under the open boundary condition. Correlation functions of both 
systems are very close to each other. We also get similar results for singlet 
pairing and spin-spin correlation functions of the same systems. 

These results show that low energy behavior of t-J ladder is approximately 
described as non-interacting system. It also indicates that the correlation 
exponent K p (— K p ) is close to unity. The direct numerical calculation of 
K p by eq.(4) shows K p ~ 1.0 for iV s = 18 system with iV e =4. It reminds 
us that the wave function of a single chain t-J model at t = J/2 is close to 
projected Fermi liquid^] and the exponent K p becomes unity in the low 



density limit |I7|. Therefore, the properties of t-J ladder are close to that of 
the single band model. It suggests that the existence of upper band E + {k) is 
irrelevant to the low energy behavior of t-J ladder in the low electron density 
region flT6 . 



5 



Next, we consider the case of high electron density. For J = 1, we can also 
compare a behavior of correlation functions and the value of K p {= 2K P ). Re- 
cently, Hayward et al. calculated correlation functions by using the density- 
matrix renormalization-group method at J/t — 1 and n = 0.8. Their results 
show that a pairing correlation is longer range than other correlations, decay- 
ing slower than ~ Rr 1 . The power of the density- density correlation function 
seems to be —2 or slightly larger. 

If we assume K p ~ 1.5 in this case, the powers of correlation functions 
can be explained with the 'duality' relation which is introduced by Nagaosa 
and Oshikawa pOf ; the powers of paring and density- density correlations are 
obtained as ~ 0.7 and as ~ 1.5 respectively. It is consist with the value of 
K p ~ 1.4, which is obtained from the Luttinger liquid relation for the N s = 20 
site system with 16 electrons. These results suggest that the Luttinger liquid 
relation holds for t-J ladder. 

Furthermore, we calculate the pairing correlations and density-density 
correlations for another system. Figs. 3 show the rung-rung pairing correla- 
tion function C SC {R) =< A A^ >, where A^ = (c^^c^,; - c^c^) 
and the density- density correlation function Ccdw(R) for N s = 16 site sys- 
tem at n=12/16. For J = 0.5, both correlation functions seem to decay as 
~ as shown in Fig.3(a). For J = 1.5, pairing correlations decay slower 
than the density- density correlations; the former seems to decay as ~ R~ ' 5 
or more rapidly and the later seems to decay as ~ R~ 2 as shown in Fig.3(b). 

Although the system size is not sufficiently large to get exponents of 
correlation functions precisely, the powers of both correlation functions seem 
to satisfy the duality relation. At least, the behavior of correlation functions 
suggests that the paring correlations dominate others for J/t>0.5. Using 

eq.(6), the critical exponent K p is obtained as ~ 1.0 and ~ 2.0 for J = 0.5 and 
J = 1.5 respectively. It consists with the behavior of correlation functions 
and indicates the validity of the Luttinger liquid relation for two-chain model. 

In the limit of J — > 0, the situation is complicated. Nagaoka ferromag- 
netic state is known to appear at the ground state of finite systems^]. How- 
ever, a very small J lifts the ferromagnetic state and changes the ground state 
to singlet. If the Nagaoka ferromagnetic phase exists in the thermodynamic 
limit, the phase may be restricted to the region with very small J. Thus, the 
lowest singlet state is not the true ground state at J = 0; nevertheless, we 
single it out as the relevant state at very small J. 
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Fig.l shows that the band structure of t-J ladder at J = is close to 
that of non-interacting spinless fermions in the ladder model. In Fig. 4, we 
compare the density-density correlation functions Ccdw(R) °f t-J ladder 
with that of spinless fermion ladder. It shows that both correlation functions 
are close to each other. Using the equations (6), the correlation exponent K p 



is estimated as ~ 0.4 |2T], which is almost consist with the value of K p (— 0.5) 



of the spinless fermion system. The same situation is already found in a 



single chain t-J model for very small J[22|. This result indicates that we 



should apply not the relation for two-chain model but single-chain model 
contrary to expectation. Probably, only the upper band E + {k) is relevant to 
the low energy behavior and charge degree of freedom is reduced to that of 
spinless fermion system in the high density region at J = 0. 

Finally, we examine the phase diagram of superconducting state in the 
J-n plane with the exponent K p . In Fig. 5, we show the possible phase 
diagram of t-J ladder. When the K p (xc) diverges, the uniform state becomes 
unstable and phase separation occurs. The boundary of phase separation 
agrees with the results of previous works [0, |9|]. Fig. 5 shows that the region 
of superconducting state appears at lower value of J ft than a single chain 
case at high electron density. It seems to resemble the behavior in two- 
dimensional case. On the other hand, the behavior of the superconducting 
region at small electron concentration resembles that of a single chain [|17] . 
Hayward et al.|14| also obtained the phase diagram in the J-n plane, but 
their value of K p is smaller than K p by a factor 2. If the value of K p is 
doubled in their phase diagram, the result of K p agrees with ours. 

The author gratefully acknowledges many useful discussions with Takami 
Tohyama. The computation in this work has been done using the facilities 
of the Supercomputer Center, Institute for Solid State Physics, University of 
Tokyo. 
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Figure captions 



Fig.l. The chemical potential /i as a function of the filling n for various 
values of J. fi is calculated for N s = 12, 14, 16, 18, and 20. The dashed 
line represents a non-interacting fermion ladder band and the broken line 
represents a spinless-fermion ladder band. 

Fig. 2. The density-density correlation functions Ccdw(R) of t-J ladder 
and non-interacting fermion systems. Ccdw(R) is calculated for the N s = 18 
sites system with 4 electrons. The dashed line has slope -2. 

Fig. 3. The rang-rang paring correlation function Csc{R) an d the density- 
density correlation function Ccdw{R) of the N s — 16 sites system with 12 
electrons (a) for J/t = 0.5 and (b) for J/t = 1.5. 

Fig. 4. The density-density correlation functions C(R)cdw of t-J ladder 
and non- interacting spinless-fermion systems. C(R)cdw is calculated for the 
N s = 18 sites system with 16 electrons. The dashed line has a slope -1.5. 

Fig. 5. A possible phase diagram of the t-J ladder as a function of J/t and 
electron density n = N e /N s . The circles represent results of N s = 16 sites 
systems. 
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